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We investigate the transport and coherence properties of a double quantum dot coupled to a single 
damped boson mode. Our numerically results reveal how the properties of the boson distribution 
can be steered by altering parameters of the electronic system such as the energy difference between 
the dots. Quadrature amplitude variances and the Wigner function are employed to illustrate how 
the state of the boson mode can be controlled by a stationary electron current through the dots. 
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I. INTRODUCTION 

The two-level system coupled to a single bosonic mode 
(Rabi Hamiltoniani) is probably one of the best studied 
models for the interaction of matter with hght^. Cavity 
quantum electrodynamics is an example where the cou- 
phng between atoms and photons can be studied in detail 
and used in order to, e.g., transfer quantum coherence 
from light to matter (control of tunneling by electromag- 
netic fields^) and vice versa^*^. 

Quantum optics usually deals with closed two-level 
systems where the total electron number on individual 
atoms remains constant and does not fluctuate during 
the interaction with the photon. This restriction can be 
lifted in semiconductor quantum dots (artificial atoms) 
by tunnel coupling to electron reservoirs. For example, in 
semiconductor cavity quantum electrodynamics^ the in- 
teraction of light with excitons can be steered by 'pump- 
ing' the cavity by resonant tunneling of electrons and 
holes. 

In this paper, we examine a single bosonic mode inter- 
acting with a system of two (bound) electronic states 
which are themselves coupled to a continuum of free 
electrons. As a concrete realization we investigate the 
stationary electron transport through a double quantum 
dot coupled to electron reservoirs and a single photon 
or phonon mode. The main idea is to control the den- 
sity matrix of the coupled dot-boson system by external 
parameters such as the reservoir chemical potentials and 
the tunnel couplings. These parameters (which in experi- 
ments can be controlled via gate voltages) then determine 
the reduced density operators of the dot and the boson 
and allow one, e.g., to modify the state of the boson by 
driving a stationary electron current through the dots. 

One of the motivations for our work are experiments 
in double quantum dots where the coherent coupling be- 
tween classical light (microwaves) and electrons can be 
detected in electronic transport i2iiSiiiii^. The open two- 
level-boson model describes the full quantum version of 
these systems with the light replaced by a photon field 
with its own quantum dynamics. 

On the other hand, quantum effects (like spontaneous 
emission) relevant for transport through double dots 
have been found to be due to phonons rather than pho- 



tons. In fact, the importance of electron-phonon cou- 
pling for transport in coupled dots is well-established by 
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Suggestions have been madi 
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explore 'semiconductor phonon cavity QED' in nanos- 
tructures where phonons become experimentally con- 
trollable. In contrast to standard semiconductor cavity 
QED^, boundary conditions for vibration modes lead to 
extremely non-linear phonon dispersion relations. For ex- 
ample, van-Hove singularities in the phonon DOS show 
up at certain frequencies^^, signatures of which seem 
to be relevant for transport through quantum dots cm- 
bedded into free-standing structures2£. Such a situa- 
tion would then be described (within a strongly idealized 
model) by the coupling of a single frequency boson mode 
to a few- level quantum dot. 

The influence of 'nano-mechanical' vibrational proper- 
ties on single electron tunneling in fact has emerged as a 
whole new area of mesoscopic transport, triggered by the 
possibility to explore electron transport through individ- 
ual moleculeaSiiS2iSi2li2^ or in freestanding2Si2L22i2S and 
movabla22i2iiSi2^ nano-structures. Phonons then are no 
longer a mere source of dissipation but become experi- 
mentally controllable with the possibility, e.g., to create 
phonon confinemenbiSiMiSi^ or the analogon of quan- 
tum optical phenomena such as coherent and squeezed 
phonon statesi^£i2Si^2iia 

The paper is organized as follows: in section m we 
introduce the model and derive a master equation for 
the density operator. We discuss transport properties in 
section Hill where a comparison to analytical solutions is 
made. In section Hvl the reduced boson density opera- 
tor, its Wigner function and the fluctuation properties 
of the boson mode under a stationary electron current 
are discussed. Finally, we conclude with a discussion in 
section M 



II. MODEL AND METHOD 

We consider an idealized situation of a single cavity 
boson mode (photon or phonon) coupled to a two-level 
electronic system (dot) which itself is connected to ex- 
ternal electron reservoirs. Our aim is to determine the 
stationary reduced density operator p{t) of the dot-boson 
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system for large times t, treating the dot and the boson 
on equal footing. The single bosonic mode (photon or 
phonon) interacts with the electrons within the dots. We 
set h=l in the following. 

A. The Hamiltonian 

An artificial 'open' two-level system can be realized 
by two quantum dots (which for simplicity we call 'left' 
and 'right') coupled to each other by a tunnel barrier 
and to a 'left' source and 'right' drain electron reservoir 
with chemical potentials /xl and ^_r. In the following, 
we are only interested in the regime of high bias voltage 
V = HL — fJ'R > We assume that the charging energy 
U required to add an additional electron to the double 
dot is much greater than V (strong Coulomb blockade 
regime). Therefore, electrons tunnel through the struc- 
ture only from the left to the right. Only one additional 
electron can tunnel into the double dot at a time, and 
the effective Hilbert space of the electronic system can 
be defined by three state3i2ii& 'empty', 'left', and 'right', 
|0) - \L) = \NL + hNR), \R) = \Nl.Nr + 1). 

The energies of the two states with an additional electron 
in the left (right) dot are denoted as El [£r)- Higher ex- 
cited states for both dots are assumed to lie outside the 
energy window 

^J■L £l,£r:> fj-R, (i) 

which defines the regime where nonlinear transport does 
not significantly depend on and ^r at low tempera- 
tures. This corresponds to the standard situationi2*i^4i& 
where only the two lowest hybridized electronic states in 
the double dot contribute to transport. 

The Hamiltonian Hd of the double dot is given in 
terms of the dot-operators til ■= \L){L\, ur := \R){R\, 
and a term Ht describing the tunneling between the left 
and the right dot, 

Hd = HA + HT^eLfiL+eRnR + TrXP + P^), 
P = \L){Rl P^ = \R){L\. (2) 

The tunnel matrix element Tc is used to describe the 
strength of the tunnehng process. 

The Hilbert space of the bosonic system is spanned by 
the usual number or Fock eigenstates |n),n = 0,l,2,3,... 
of the harmonic oscillator 

Hp = Lua^a (3) 

with frequency u and creation/annihilation operators a 
and fulfilling [a, a^] = 1. The coupling between the 
electronic and the bosonic system is described by four 
microscopic coupling constants. 

Hep = (miL + f372R + -f*P + jP^){a + a^). (4) 

Here, we assume the coupling constants a, (3, and 7 as 
given parameters. Their precise form can be calculated 



from microscopic details such as the many-body wave 
functions of the dot electrons. It should be noted that the 
interaction between the bosonic system and the reservoirs 
is not considered here. 

The two electron reservoirs are described by 

Hres = ^kclck + ^kAdk, (5) 

A: k 

where the sum is over all wave vectors k in both the 
left (L) and right (R) reservoirs and spin polarization is 
assumed for simplicity . The coupling between dot and 
reservoirs is described with two dot operators 

SL = |0)(L|, SR^\0){R\ (6) 

and the tunnel Hamiltonian 

Hv^J2^Vi^^lsL + V,''dlsR + c.c) (7) 

k 

with tunnel matrix elements Vjf and Vj^. 
The total Hamiltonian is written as 

H = Hq + Ht + Hy + Hep, (8) 

where Hq = Ha + Hp + Hres is the free Hamiltonian in 
the interaction picture introduced below. 

B. Equations of Motion 

The density matrix of the total system (dot, boson 
and reservoirs) is given by the Liouville-von-Neumann 
equation and formulated in the interaction picture with 
respect to H^. For the non- linear transport window Eq. 

considered here, the chemical potentials in the equi- 
librium reservoirs are such that the Fermi distributions 
in the left and the right lead can well be approximated by 
/^(e) = 1 and /^(e) = 0. We treat the reservoir couphng 
term, Hy, in second order perturbation theory, neglect- 
ing Kondo physics^ throughout so that transport can be 
described by a master equationiSiiSilSii^. Transforming 
back into the Schrodinger picture produces the following 
master equation for the reduced density operator of the 
system (dot -I- boson), 

j^pit) = -l[HA+Hp + HT+Hep,p{t)] (9) 

- ^{sLslp{t)^2s{p{t)sL+p{t)sLsl) 

- —(sRSRAt) ~ '2-SRp{t)s^jj + p{t)s^j^SR) 

- '^{2apa' — a) ap — pa) a) 

The perturbation from the reservoir is given by the tunnel 
rates 

T„R ^ 2.Y\yk"'m^L/R - el'\ (10) 

k 
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FIG. 1: Stationary current through double quantum dot cou- 
pled to single boson mode with Tl ~ = 0.1, Tc = 0.01, 
7b — 0.05 and varying 2g = a = —[3. (electron-boson cou- 
pling, cf. Eq. l(4j). A'': number of boson states. 



In the last line in Eq. (Q, we introduced a term that 
describes damping of the bosonic system at a rate 
corresponding to photon or phonon cavity losses. 

Taking matrix elements of the master equation, Eq. 
, one obtains a system of linear equations for the ma- 
trix elements of the density operator which are given ex- 
plicitly in appendix ^ In the stationary state of these 
equations, the time derivatives are zero, and a coupled 
set of linear equations is found. In order to numerically 
solve these equations, the bosonic Hilbert space has to 
be truncated at a finite number N of boson states. The 
total number of equations then is 5A^^-|-10iV-|-5, remem- 
bering that there is always an equation for n = m = Q. 
The numerical solution becomes a standard matrix in- 
version, though book-keeping of the matrix elements has 
to be done carefully. The trace over electron and bo- 
son degrees of freedom has to to obey the normalization 
condition 

Y,[{nAPstat\n,Q) + {n,R\pstat\n,R) (11) 

n 

+ {n,L\pstat\n,L)] = 1. 

We mention that without this condition, the matrix be- 
comes singular, and cannot be solved. 



III. TRANSPORT PROPERTIES 

The first quantity of interest that one can obtain from 
the matrix elements is the stationary electron current 
through the double dot. 



A. Stationary Current: Numerical Results 

In the stationary state, the current operator is defined 
via the rate of electrons tunneling from the left dot to 
the right dot, 

i = Q-nL = iTc{P^ -P) + i{-fP^ - j*P){a + at). (12) 

Here and in the following we set the electron charge — e = 
1. For a boson mode with wave vector Q, the electron- 
boson coupling constants can be expressed in terms of 
matrix elements of the left and right dot states, a = 
A(L|e^Q'-|L), (3 = XiRle^^^lR), 7 = A(L|e*'^''|i?), where 
A is the microscopic constant for the interaction of the 
boson with the electron in 2D. 

Our formalism works for arbitrary choices of coupling 
parameters a, (3, 7, but for simplicity we restrict our- 
selves to simple, non-trivial cases for the numerical cal- 
culations. We set 7 = which corresponds to neglect- 
ing non-diagonal terms that for relatively sharply peaked 
electron densities in the dots only weakly contribute^. 
For a sharp electron density profile, one furthermore has 
(3 = ae**^'* where d is the vector connecting the left and 
right dot centers. Identical energy shifts in both dots, 
corresponding to /3 = a, have no effect on transport. In 
the following, we choose Qd = tt and therefore electron- 
boson coupling constants 

a = -/?=2.g, 7-0. (13) 

Furthermore, the tunneling rates between the reservoirs 
and the dots are generally kept smaller than the ener- 
gies of the dots and boson states. This produces clear 
and sharply defined transport characteristics. We use 
the boson frequency as an energy scale and set uj = 1 
in the following. Using the typical phonon cavity energy 
huj w 100 /ieV from a recent experiment with a freestand- 
ing phonon cavitjs^S yields a typical electron current scale 
of eoj = 24 nA in Figs. [l]and||| 

Fig. n illustrates the stationary current as a function 
of the energy difference e. For zero boson coupling g = 0, 
we reproduce the Lorentzian resonant tunneling profile 

Istat = T^Tr/[T^{2 + Tr/Tl) + Tl/A + £2] (14) 

as first derived by Stoof and Nazaro\ii£. 

For finite electron-boson coupling and positive e, reso- 
nances appear when 

e = El - sr^ noj, (15) 

i.e., the energy gap of the dots becomes equal to multiple 
integers of the boson energy. At these resonances, elec- 
trons tunneling from the left to the right dot can excite 
the boson system. Note that for larger Tc, e = el — £r 
in Eq. (|15|) has to be replaced by the energy difference 
A = \J e'^ + 4T2 between bonding and anti-bonding state 
in the double dot. 
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N=10,Yb=0.05 
POLYb=0.05 
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FIG. 2: Current as in Fig. for coupling 2g = 0.1 and 
different boson damping 71,, alongside corresponding polaron 
transformation result, Eq. 1191 



The current profiles with differing damping rates, 
Fig.©, illustrate that the system is extremely sensitive 
even to very small 'damping'. For small but finite bo- 
son damping, we found that two profiles for N = 5 and 

= 10 match extremely closely. This illustrates that a 
finite boson damping removes numerical problems due to 
the truncation of the boson Hilbert space at a finite N. 
As expected, there are no peaks at all on the absorption 
side of the profile (e < 0): the damped boson relaxes to 
its ground state. 



B. Perturbation Theory in Tci Polaron 
Transformation Method 



An analytical expression for the stationary current 
can be obtained from a polaron transformation and a 
polaron-transformed master equation that is valid for 
small coupling Tc between the two dotsi^. After tracing 
out the bosonic degrees of freedom one obtains an equa- 
tion of motion for the time-dependent expectation values 
of the operators nj^/ji and P'^^ which can be transformed 
into Laplace-space, 



T, 



{l/z-(ni)(z)-K.)(z)} 



T, 



r. 



^^{{P)[z)-{P^)[z)]-^{n!,){z) 
-iT, {{nL){z)Ce{z) ~ {nR){z)C%{z)] 



-{P){z)C,{z). 



(16) 



Here, the coupling to the boson system enters through 
the correlation function 



Ce{z) 



dte 



-zt-\-iet 



C{t) 



C{t) = {X{t)X^{t^O)), X^D 
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(17) 



of the displacement operator of the bosonic mode, 

l?(0=cxp(eat-ra). (18) 

The equations 1)16(1 can be solved algebraically. One then 
obtains the expectation value {I)t-^oo of the stationary 
current from the 1/z-coefficient of the /(z)-expansion 
into a Laurent series^ for z — s- 0, 



{I)t 
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(19) 



where Ce = C,{z^Q). 

The correlation function, C{t), enters via a factor- 
ization assumption in the polaron-transformed master 
equation and has to be calculated from an effective den- 
sity operator of the bosonic mode. For dissipative 
pseudospin-boson systems with a continuous spectrum 
of modes and a thermal equilibrium for pe, reasonable 
results can be obtained within this method for small 
coupling, in particular in comparison with perturbation 
theory^^*^. However, in the present case with just one 
single mode it is not clear whether or not a factorization 
into an effective dot and an effective boson density opera- 
tor yields reasonable results, and in principle one should 
at least determine pb self-consistently. For small cou- 
plings g, however, the assumption of a decoupled time- 
evolution of the boson should work well, and we proceed 
with a master equation for pe- 

d 



dt 



PB{t) = -i[uja^a,p^] 



(2apBa 



a apB ^ Pbci ci 



(20) 



that describes the free time-evolution of the boson in 
presence of boson damping with damping constant jf, > 
0. Assuming the boson initial condition to be its ground 
state, one obtains 



Cit)=exp[-\e{ 



}, e = ^ 

J UJ 



(21) 



Fig. (jSJ shows that the analytical expression for the cur- 
rent, Ea. (|19|l . compares quite well with the numerical 
results. In particular, this demonstrates that the formu- 
lation in terms of a general boson correlation function 
C{t) works reasonably well not only for boson systems 
with continous mode spectrumiSii^ but also for the sin- 
gle mode case discussed here, at least for small coupling 
constants g. We point out that the existence of a finite 
boson damping is crucial here and that for larger g, the 
comparison becomes worse, as expected. 
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IV. BOSON DISTRIBUTION 

Turning now to properties of the bosonic system, our 
primary interest is how the boson mode can be controUcd 
via parameters of the electron sub-system such as the bias 
£ and the tunnel couplings Tyi/l- In the stationary state, 
the reduced density operator of the boson is 



Pb = lim Trdotp(0 = hm p"{t). 



(22) 



1=0, L,i? 



Boson state detection 



Before discussing the boson state pb we describe a pos- 
sible experimental scheme to directly detect its proper- 
ties. The basic idea is to use another nearby double quan- 
tum dot as detector of the boson, similar to the scheme 
for detecting quantum noise in mesoscopic conductors as 
proposed by Aguado and Kouwenhoven^. 

The detector consists of a double quantum dot very 
similar to the one discussed above. The boson state gives 
a contribution to the inelastic, stationary current 







Gd{t) 



Tr 



PbXd{t)X\ 



(23) 



through the detector in lowest order in the detector tun- 
nel coupling Td. Here, E denotes the detector dot energy 
difference, D again denotes the unitary displacement op- 
erator, Eq. H18|l . gd the coupling constant in the detector 
double dot, and ph the boson state. Eq. describes 
the detector current I{E) resulting from a given boson 
state Ph as calculated in the next section. 

Before discussing our numerical results for pb, it is illus- 
trative to look at a few special cases. For example, if the 
boson state p^ describes a thermal equilibrium at tem- 
perature T, one obtains the inelastic current spectrum 
with absorption E < Q and emission E > Q branchesi^, 
as has been observed by Fujisawa et alii^ for the case of 
equilibrium (multimode) phonons. 

If Ph = \n){n\ describes a pure, n-boson number state, 
the function Cd{t) in the detector current, Eq. (|23|l . be- 
comes {£,= ^gd/uj) 

Cd{t) = (n|I?(ee-*)i?(-C)|n) - 



(24) 
|2^ 



where is a Laguerre polynomial, cf. appendix^ Here 
we assume that the boson time-evolution is undamped 
in the detector, i.e. governed by the boson Hamiltonian 
Hh = uia^a. Expanding up to lowest order in ^ = 4:gd/uj 



and using L'!^{x) = 1 
Eq.lESll and Ea. (|^ 



0{x ), one obtains from 



Id{E) « T!2n [{1 - \e{2n + 1)} <5(i?) (25) 
+ \^\^{n5{E + ij) + {n + l)S{E-uj)}] 




n=0 
n=2 
n=4 
n=6 
n=8 
Norm 
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e=eL-eR (100|xeV) 

FIG. 3: Boson profiles n=0,2,4...10 for jb = 0.0001, iV = 10, 
2,g = a = = 0.3, 10 x Tc = Fl = Fa = 0.1. 



Clearly, the height of the inelastic current peaks is de- 
termined by the quantum number n, i.e. the absorption 
peak at = —uj < scales with n, and the (stimulated) 
emission peak at i? = a; > scales with n + 1. 

The above example can be even generalised. In ap- 
pendix El we show how to reconstruct an arbitrary bo- 
son state Ph = J2nm Pnm\n) {m\ from the function P{E) 
(or equivalently the stationary detector current spectrum 
Id{E), Eq. H23|l 'l. This demonstrates that the properties 
of the boson state discussed in the following can be di- 
rectly related to an experimentally accessible quantity. 



B. Numerical results 

We first analyse the boson mode via moments of its dis- 
tribution function such as occupation number and Fano 
factor, before discussing the full distribution function in 
the Wigner function representation. 

The boson occupation probability p„ = {pb)nn is the 
probability for the boson to be in one particular number 
state n. Furthermore, the boson Fano factor F'^ is calcu- 
lated from the boson occupancy as the variance of boson 
number over the mean number of bosons. 



Similarly, the variance Aa;^ = (a:^) — (x) 
(p^) — (p)^ of the position and a momentum coordinate 



(26) 



and Ap"^ = 



X = T= — , p = 



x/2 



^/2 



(27) 



can be obtained from ph, cf. Eq. ljA7|l in appendix 1X1 
They provide an indication of the quantum fluctuations 
in the boson states, and prove useful in comparison with 
the Fano factor profile. It should also be apparent that 
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-2-10 1 2 
e=EL-eR (100 ^eV) 

FIG. 4: Fano versus s profile for 2g = a = -/3 ^ 0.1, TV = 10, 
7b = 0.05, 10 X Tc = Fl = Ffl = 0.1. Inset: Comparison be- 
tween coiierent state, Eg. 12911 . and numerical result for boson 
distribution pn at e = —4. 



We have checked that the probability p„ fulfills the 
normaUsation J2nPn = 1- Foi' accurate results the oc- 
cupation probability should tend to as n — > A'^, where 
N is the dimension of the truncated boson Hilbert space. 
This condition is easily achieved in presence of a finite 
boson damping rate, jb > 0. In this case, it is possible 
to obtain a boson distribution which is centered around 
the lower boson number states, and is excited when a 
resonant interaction with the electron occurs. 

It is more illustrative to first discuss a weakly damped 
case, as shown in Fig. |3 We can see quite clearly how at 
the positive resonance energies, the occupation probabil- 
ities of the states spreads into the higher number states. 

The Fano factor F'^ illustrates fluctuation properties of 
the boson mode with F'^ = 1 corresponding to a Poisso- 
nian boson number distribution and F'^ < 1 {F^ > 1) to 
a sub (super)-Poissonian distribution. Fig. 0] illustrates 
the Fano factor vs. e for the strongly damped case. 

This profile suggests the reduced boson states are co- 
herent (Poissonian) for e < 0, and super-Poissonian at 
the resonance energies, e ^ lu,2lu, .... 



C. Coherent state and x-p variances 

For £ <C 0, nearly no current is fiowing, the electron 
is prcdominatly localized in the left dot, and one can 
approximate the operator az = \L){L\ — \R){R\ by its 
expectation value (az) = 1. Then, the boson system is 
effectively described by 



(28) 



which is a shifted oscillator, the ground state \GS) of 
which is the coherent state \GS) = | — 2g/(jj) with a\z) = 
z\z). This can be seen by introducing new operators 



Position Variance 
Momentum Variance 
(Ax.Ap) 




0.5 1 1.5 
8=eL-eR (100 (teV) 
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FIG. 5: Position variance, momentum variance, and an esti- 
mation of the combined uncertainty. 2g = a — —(3 — 0.1, 
N = 10, 76 = 0.05, 10 X Tc = Fz, = Ffl = 0.1. 



b := a+2g/u! whence Hcs 
In this regime, we have 



iub'fb-Ag^/uj and6|G5') = 0. 



Pb^\z){z\, z^-2g/u, 



(29) 



and p„ 



|(n|G5)p is given by a Poisson distribution, 
Pn — |2;p"e~l^l /nl. One can plot this against the nu- 
merically obtained p„ and check that the distribution is 
indeed Poissonian for small £ <C 0. This is represented 
in the inset of Fig. 0]in detail for the first 4 boson Fock 
(number) states. 

The quantum fiuctuations of the boson system is de- 
scribed by the variance of the position and momentum 
operators. Fig. |S1 presents a direct comparison of the 
position and momentum variance of a strongly damped 
system. The uncertainty principle holds, AxAp > i, and 
thus Ax^Ap^ > i. The state with minimum uncertainty 
occurs as expected, for el < e_R, with Ax^Ap^ = j. This 
is a reassuring result, as it suggests both that the boson 
is coherent in this area, and that the Fano factor's Pois- 
sonian distribution for this coherent state is correct. We 
mention that since both Aa; > 0.5 and Ap > 0.5, we 
find no squeezing of the boson mode for the parameters 
checked here. 



D. Wigner Function 

The Wigner function of the density operator pf, for a 
bosonic mode is a representation of pi, in x-p-spa.ce, 
cf. Eq.lEII). It is defined as^ 



W{x,p) -Tr {pbD{2a)Uo) 

TT 



ip 



V2 



(30) 



where D{a) :~ exp[aa^ — a* a], Ea. H18|) . is a unitary dis- 
placement operator and Uq :~ exp[i7ra^a] parity opera- 
tor for the bosonS. W{x,p) is a symmetric Gaussian for 
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£ =0.00 



= 0.25 



= 0.50 



DISCUSSION 




FIG. 6: Wigner distribution functions for the bosonic mode. 
Parameters are Vl = = Tc = 0.1, 7;, = 0.005, g — 
0.2, A'^ = 20. Stationary current resonances occur at £ = 
0.0,1.0,2.0,.... 



a pure coherent boson state and a symmetric Gaussian 
multiphed with a polynomial for a pure number stated. 
In our model, the shape of W{x,p) therefore indicates 
how close to these limiting cases the actual stationary 
state pb of the boson, Ea. lf^ . is. In particular, this is 
a convenient way to represent the 'steering' of the bo- 
son mode when external parameters (like £ or T n/i^) are 
changed. 

Using the Fock state basis {|n),n = 0,1,2,...} and 

'^o = Er=o(-l)"l")H, wefind 

^ oo 

W{x,p) = - {~ir{n\pb\m){m\D{2a)\n). (31) 

n, 771—0 

It is useful to split the sum into diagonal and non- 
diagonal parts, to use (TO|pf,|?i) = (?i|pf,|m)* and [m > n) 



(TO|£'(Q;)|n) = \ — -c 



'5l«l'L™-"(|«|2) (32) 



= {m\D^a)\72)* = (-l)'"-"(m|D(a)|7i)*, 

is a Laguerre 



where again a = {x + ip) / \/2 and L™ " 
polynomial, cf. appendix IbI This leads to 



W{x,p) = -f] 



oo 



oo 



E E 

n— m— n+1 



{-\T(n\Ph\n)l,n\D{2ot)\n) (33) 
(-l)"2Re[(?i|p6|?7i)(m|i:)(2a)|n)] . 



Using our numerical results for the stationary state 
p of the coupled system, we easily obtain Wigner func- 
tions from Ea. (|22|l and Ea. H33() for different values of 
the energy difference e between the two dots as shown 
in Fig. Between two resonance energies e = nuj, 

W{x,p) closely resembles a Gaussian. At and close to 
the resonance energies, the distribution spreads out in 
rings around the origin, which is consistent with the in- 
creased Fock state occupation numbers and the increased 
position and momentum variances at these energies. 



Our investigation is based on a numerical method 
which allows the density matrix of the boson-dot system 
to be solved for arbitrary dot parameters (e, Tc, r^/zj) 
and boson damping rate 7^. Since we have to truncate 
the boson Hilbert space, the electron-boson coupling con- 
stant g has to be restricted to small values 5 < 1. In 
contrast to the determination of the (pure) eigenstates 
of an isolated dot-boson system (Rabi-Hamiltonian) , the 
numerical effort for our mixed state (density operator) 
is much bigger here. Although not discussed in this pa- 
per, we suggest that the strong coupling regime could be 
reached numerically by a polaron transformation of the 
master equationi^ without the factorization assumption 
employed in section IIII Bl Fortunately, at present the 
small coupling regime seems to be valid for experimental 
situations with quantum dot9iiiiiiSii&. 

Our results suggest that there is no resonant interac- 
tion on the e < side of the current peak as long as 
the boson system is damped and any many-body excited 
electron states can be ignored. For e > we have found 
strong excitations of the boson mode occuring at reso- 
nances given by multiples nuj of the boson frequency uj. 
These correspond to the photo-satellite peaks in reso- 
nant tunneling in electromagnetic AC fields^iiS which in 
the classical case, however, appear on both sides of the 
main resonance n = in contrast to the quantum case 
considered here. The resonances at e = nu can thus be 
interpreted as the emission of phonons (photons) by the 
electron as it tunnels through the dot. It should be pos- 
sible to detect them in the stationary current through 
double quantum dots in 'phonon resonators'. With in- 
creasing electron-boson coupling, the visibility of these 
side peaks increases whereas the main resonance at e = 
is reduced. 

We have discussed that an additional double dot can 
serve as a detector of the stationary boson state pb, if 
its energy difference E is tuned independently. In the 
single mode case considered here, the stationary detector 
current spectrum Id{E) is uniquely related to the com- 
ponents of ph in the number state basis. In order to 
characterize the stationary boson state itself, we have in- 
vestigated the Fano factor and the quantum fluctuations 
in the quadrature amplitudes of the boson. For e <C 
the boson state is perfectly described as a coherent state 
whereas for e ^ the interaction with the quantum dot 
produces occupied excited states in the boson mode that 
could best be visualized by a broadening of the phase 
space (Wigner) distribution function of the boson. 

Although in this paper we have only calculated station- 
ary properties, our method can be extended to calculate 
the noise-spectrum via the quantum regression theorem 
within the master equation framework as well^. It re- 
mains a task for the future to analyze the relation be- 
tween the boson quantum fluctuations and the current 
noise in detail. 
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APPENDIX A: MATRIX ELEMENTS OF THE 
DENSITY OPERATOR 



The matrix elements of the dot-boson density operator 
p are defined as 

Pnm '■= {n,i\p\i,m), i = 0,L,R 

Pnm ■= {n,i\p\.hm), i,j = Q,L,R, (Af) 



IS 

d 



^Pnm i[Pnra[^irn - n) + {e L - £/?)] 



(A5) 



+ Vm + Ip^^+i + V^Pn^-l) 



nm+1 
BR LL 



■7n 



/mp„ 

Pnm/i t2 [Pnmi 



+ Y^'2Pn+l,m+lV^^Vm + 1 - + m)]. 

Taking the trace over all the boson and electron states, 
the expression Eq. ((T^ for the electron current operator 
reads 



71 



(A6) 



where n and m refer to boson Pock states. 

The matrix elements for the empty dot state obey 



dt' 



[t^{m - n)plJ - Tl [p?„„] + TrIp^J (A2) 



+ y[2p° +!,,„+! + iVm + 1 - p"„,(n + m)], 



those for the right dot state 



dt' 



^ i[u:{m-n)p^^+T,{p^^^-p^^) (A3) 

+ /3(V^Pmn-l + Vm + IPnm+i 

- V^Pn-lm - V?l + IPn+lJ) 

+ 7*(V^P,™-1 + \/'«+lpnm+l) 



+ y[2p^+i,m+iVnTTV^^rfT-p^„(ri + m) 



- 7* {Pnn-l + Pni+1 \/^)] • 

The variances of the boson position and momentum co- 
ordinate, Ea. H27() . are obtained by performing the trace 
over the dot variables 0,L,R, 



Ax^ = ^Trdot y^(%/nV^ 



n - lpnM-2 



(A7) 



Ap2 



(2n + l)pn,n + ^/n+ Wn + 2p„^„+2) 

n,n+lJJ ; 

n 

^Trdot y^(-VnV?i - lp„,n-2 



(A8) 



+ {2n + l)pn,n - \/n + lVn + 2/9„^„+2) 

- -[Trdot y^}~iVnpn,n-i + Wn + lp„,„+i]]^. 



APPENDIX B: DISPLACEMENT OPERATOR, 
COHERENT STATES 



and for the left dot state 



dr 



+ a{-VnTTpn+lm - VnPn-lm 



(A4) 



+ 7(\/^P™^-l + V^^Pm) 



Tl [Pnr 



-[2a 



n+l,rn+ 



iV"- + iVm + 1 - p„„j (n + ; 



The equation of motion for the off-diagonal elements 



Here, we summarize some useful properties of the uni- 
tary displacement operator 



D{z) = e^" 

— 1 I 

^ e 2i 



{dHz))-'=D^{-z) 



where z is a complex number and we used the operator 
exponential e^+^ = e^e-^e-(i/2)[A.B] fo^. [[A,B],A] = 
[[A, = 0, cf. also^. A coherent boson state \z) is 

defined as eigenstate of the annihilation operator. a\z) ~ 
z\z), where z is a complex number. It can be generated 
from the boson vacuum |0) as 



\z)^D{z)\0). 



(B2) 
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Extremely useful is the relation 



(B3) 



for arbitrary complex numbers a, (3. Coherent state ma- 
trix elements of D{a) follow as 



{l3\D{a)\P) = e-5l"l'e2»s(ar). 



(B4) 



Number state matrix elements can be obtained using 
\n) EE (1/V^)(at)"|0), 



{m\D{a)\n) = e^'"'' x 



(B5) 



|g(2i-Q*)ag(22+Q)a'^ I 



21=22=0 



With (0|e(-'i~"*)'^e(^2+a)at|Q^ ^ ^{z^~a')(^2+a) ^ for m > 
n the differentiation yields 



{m\D{a)\n) 



\J n\m\ 



.(B6) 



22+a 



Comparison with the generating function of the Laguerre 
polynomials^ yields Ea. (|32|l . A corresponding expres- 
sion can be derived for n > m. 



APPENDIX C: RELATION BETWEEN P{E) AND 
STATIONARY BOSON STATE pb 

Again we assume that the boson time-evolution is un- 
damped in the detector, i.e. governed by the boson 
Hamiltonian Hi, = oja^a. Expressing an arbitrary boson 
state pb in the number state basis, pt ~ X^nm Pnm\n){m\, 



Fourier transforming the function P{E), Ea. H23|l . yields 



Pit) EE — I dEe~^^'P{E) = Ca{t) + 



(Cl) 



Cn,n{t) = {m\Di^e'^')D{-0\n) + {m\D{OD{-^e-'^')\n). 

The functions c„,„(i) can be calculated analytically, us- 
ing Eq. I)B4|I and the matrix elements {ni\D{z)\n) of the 
unitary displacement operator, Eq. H32|l . They are pe- 
riodic in 27r/a;, and from integrating Ea. ljCnfl over one 
period one obtains a linear relation between the Fourier 
coefficients P'^ and cj;^ , 



P' 



(C2) 



pk — 



27r 



2tt 



P{t)e''"^'dt 



^nrn(t)^ ^t 



The P'^ have to be determined by numerical integration 
from the (experimentally given) P{E). Regarding (nm) 
as a single index, Ea. ljC2() is a linear equation that can 
then be solved for the coefhcients pnm of the boson state 
by inverting the matrix cjj^, the coefficients of which 
are given as Fourier coefficients of known expressions. In 
practical terms, the number of boson states as well as 
the number of Fourier coefhcients have to be restricted 
in order to make this inversion feasable. 

Eq. ljC2|l establishes the relation between P[E) (or, 
via Eq. the detector current Id ~ TjP{E) in lowest 

order Td) and an arbitrary single mode boson state ph- 
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